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Abstract. Geometrical objects with integral side lengths have fascinated mathematicians through the ages. 
We call asetP={pi,...,Pn}CZ^a maximal integral point set over if all pairwise distances are integral 
and every additional point Pn+1 destroys this property. Here we consider such sets for a given cardinality 
and with minimum possible diameter. We determine some exact values via exhaustive search and give several 
constructions for arbitrary cardinalities. Since we cannot guarantee the maximality in these cases we describe an 
algorithm to prove or disprove the maximality of a given integral point set. We additionally consider restrictions 
as no three points on a line and no four points on a circle. 



1. Introduction 

Geometrical objects with integral side lengths have fascinated mathematicians through the ages. A very 
early example is the Pythagorean triangle with side lengths 3, 4, and 5. A universal framework for most 
of these objects are integral point sets. By an integral point set we understand a set of n points in an 
m dimensional Euclidean vector space E™^, where the pairwise distances between the points are integral. 
Those integral point sets were studied by many authors, see |9 1 for an overview. From a combinatorial point 
of view for a given cardinality n and a given dimension m the question on the minimum possible diameter 
d(n, m), this is the largest distance between any two points, arises, see lfT6l [1911201 for an overview. 

To obtain some interesting discrete structures one could also require some additional properties. One 
possibility is to request, that besides the distances also the coordinates must be integral. Another classical 
possibility is to forbid subsets of three points on a line or four points on a circle. The question of P. Erdos 
whether there exists a set of seven points in the plane with no three points on a line, no four points on a circle, 
and pairwise integral distances, has recently been answered positively, see tl4l . If all three mentioned 
additional properties are required simultaneously one speaks of TLm-clusters, see ll22ll . In this article we 
request that besides the distances also the coordinates of the point sets are integral and restrict ourselves to 
dimension 2. Additionally we consider the cases where no three points are on a line or no four points are on 
a circle. 

In finite geometry one is sometimes interested in point configurations which are maximal with respect to 
some property. This means that it is not possible to add a point without destroying the requested property. 
Here we consider integral point sets which are maximal, meaning that there does not exist an additional 
point X with integral distances to the other points of the point set. 

1.1. Related work. There have been extensive studies on integral point sets in Euclidean spaces. Some 
authors also consider other spaces, e. g. Banach spaces f6l, integral point sets over rings |13|, or integral 
point sets over finite fields 12] [TTI [15] . In O the authors consider integral point sets over 71? and conjecture 
some examples to be maximal. As an answer to their open problems in |[T2l the authors describe an algorithm 
to prove the maximality of a given integral point set and prove the conjectures of 1^3]. 
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1.2. Our contribution. In this paper we describe algorithms to efficiently test integral point sets for max- 
imality and to determine possible extension points. To deal with the isomorphism problem we describe an 
algorithm which transforms a given plane integral point set into a normal form in O (n^) time, where n is 
the cardinality of the point set. We give several constructions of integral point sets over l? which have a 
given cardinality and fulfill additional conditions such as that there are "no three points on a line" or "no four 
points on a circle". Although we cannot prove the maximality of the point sets obtained with the proposed 
constructions in general, we conjecture this property for many of our constructions. By exhaustive search 
we have determined some exact minimum diameters of integral point sets over l? with given cardinality 
and with or without additional conditions. We give constructive upper bounds in most cases and conjecture 
them to be the exact values. 

1.3. Outline of the paper. In Section[3]we state the basic definitions and in Section|2]we describe the basic 
algorithms to deal with maximal integral point sets over 1? . These include an algorithm to exhaustively 
generate Heronian triangles up to isomorphism, an algorithm to determine all possible embeddings of an 
Heronian triangle on the integer grid 1? , and an algorithm that determines all points of T?- which have 
integral distances to three given points in T?- with pairwise integral distances. The last mentioned algorithm 
enables us to algorithmically prove or disprove the maximality of a given integral point set. Since we intend 
to consider integral point sets up to isomorphism we introduce normal forms of integral point sets and 
algorithms to obtain them in Section |4] We deal with the keyquestion of maximal integral point sets over 
1? with given cardinality and minimum diameter in Section p] Several constructions for maximal integral 
point sets, where the maximality is not guaranteed but very likely, are described in Section [6] In Section [7] 
we deal with additional properties as "no three points on a line" and "no four points on a circle". We finish 
with a short conclusion and an outlook in Section[8] 

2. Basics 

Definition 2.1. An integral point set over is a non-collinear set T of n points in the integer grid Z^, 
where the points have pairwise integral distances. 

For brevity we only speak of integral point sets and assume that the coordinates of the points are integral 
numbers, too. 

Definition 2.2. We call an integral point set y over maximal if for every x G Z^\CP the point set CP U {x} 
is not an integral point set. 

The existence of maximal integral point sets in the plane is guaranteed by a famous theorem of N.H.. An- 
ning and P. Erdos, respectively its proof. 

Tlieorem 2.3. An infinite set T of points in the Euclidean space with pairwise integral distances is 
situated on a line. ||Tl|4l 

Proof. We only prove the statement for dimension ra = 2, as the generalization is obvious. If A, B, and 
C are three points not on a line, we set k — max {AC, BC} and consider points P such that |PA — PC| 
and |PB — PC| are integral. Due to the triangle inequalities the attained values are in {0, 1 , ... , k}. Thus the 
point P lies on the intersection of two distinct hyperbolas, where we have at most k + 1 choices for each 
hyperbola. Thus there are at most 4(k + 1 )^ possible locations for the point P. □ 

This proof can clearly be converted into a constructive algorithm. Given three points A — (xa,iJa), 
B = ('^BiUb). and C — (xcUc) in J" C Z^, which are not on a line, the problem of determing points 
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P — (xp , ijp ) at integral distance to A, B, and C is reduced to the problem of solving the equation system 

V (xa -xp]^ + (ua -UpF - V(xc -xp)^ + (ijc -Dp)^ = di 
^(xb -xp)2 + (ye -^vY - V(xc -xp)2 + (yc -^vY = da 

where di G {—AC, . . . , AC} C Z and di S {— BC, . . . , BC} C Z. If there exists no integral solution in 
Z^XT, then the point set T is maximal. This algorithm was already used in lfT2l to prove the maximality of 
the integral point sets of Figure [T] 




Figure 1 . Examples of maximal integral point sets. 



Since this algorithm is essential for our article we will go into the details how to solve equation system 
[T| To get rid of some of the square roots we add \/(xc — xp)^ + (yc "Up)^ on both sides and square the 
expressions afterwards: 

(xa -xp)^ + (yA --yp)^ = df +2di v/(xc -xp]^ + (yc -yp)^ + (xc -xp)^ + (yc -yp)^ 
(xb -xp)2 + (yB --yp)^ = da +2d2V(xc -^vY + (yc -^vY + (xc -xp)^ + (yc -yp)"^ 
Rearranging yields 

(x^ +y^ -x^ -y^ - df) +2(xc -xa)xp +2(yc -yA)yp = 2di a/(xc -xp)2 + (yc -yp)^ 
{x| +yi -x^ -y^ - d|) +2(xc -xb)xp +2(yc -yB)yp = 2d2V(xc -xp)^ + (yc -yp)^ 

(2) 

If di =0 then the first equation corresponds to a linear equation 

cixp + C2yp +C3 = 0, (3) 

where not both Ci and C2 are equal to zero, since A 7^ C. If we square the second equation of (j2| we 
can substitute one variable using equation ([3]l and obtain a quadratic equation in one variable, which can be 
easily solved. The case, where d2 = is similar Here we use the second equation of Q to obtain equation 
(j3]l (we have Ci 7^ or 02 7^ due to B 7^ C), and substitute it into the squared version of the first equation 
to obtain the quadratic equation in one variable. In the remaining case we have di , d2 7^ 0. Here we subtract 
di times the second equation of (j2]i from d2 times the first equation of (|2| to obtain equation (|3]l (we have 
Ci 7^ or C2 7^ since the points A, B, and C are not located on a line). Now we can square one of the two 
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equations of Q and subsitute one variables using equation Q. Again we end up with a quadratic equation 
in one variable. At the end we have to check if the obtained values (xp,yp) are solutions of the original 
equation system ([TJ. 

Definition 2.4. For an integral point set T its diameter diam(CP) is given by the largest distance between a 
pair of its points. 

We remark that the left integral point set of Figure [T| has diameter 8 and the right integral point set of 
Figure [T]has diameter 25. 

3. Exhaustive generation of maximal integral point sets 

To obtain interresting examples of maximal integral point sets we utilize computers to exhaustively generate 
maximal integral point sets. In the following we will describe the algorithm used. For a given diameter d we 
loop over all non-isomorphic Heronian triangles (having integral side lengths and integral area) A — (a, b, c) 
with diameter d = maxja, b, c}. Utilizing the Heron formula 



^ _ V(a + b + c)(a + b-c)(a-b + c)(-a + b + c) 
~ 4 
for the area of a triangle we can generate this list e.g. by the following short algorithm: 

Algoritlim 3.1. {Generation of Heronian triangles) 
input: diameter d 

output: complete list of Heronian triangles with diameter d up to isomorphism 
begin 
Q = d 

for b = I'^l , . . . , a do 

for c — g + 1 — b, . ■ . , b do 

if — ^ e Z men 

output (a, b, c) 

end 

For a more sophisticated and efficient algorithm we refer to ifTSl . The next step is to embed a given 
Heronian triangle A = (a, b, c) in the plane integer grid Z^. Here we can utilize two conjectures, which are 
theorems for dimension m = 2, see e.g. fS). 

Conjecture 3.2. Let y C Q™' be a finite set of points such that the distances between any two points ofT 
are integers. In this case one can find an Euclidean motion T such that T( C 7"^. 

Conjecture 3.3. Let T C Z"^ be a finite set of points such that the distances between any two points of 7 
are integers and divisible by an integer k. In this case one can find a set J" C Z"^ such that J" • k (the set 
y scaled by a factor kj is congruent to 7. 



Since Conjecture 3.2 is a well known theorem for dimension m — 2, see e.g. Q, for every Hero- 
nian triangle A(Q,b,c) there exists an embedding in the plane integer grid Z^. We remark that there 
may be several embeddings for the same triangle A — (Q,b,c), which lead to different results. If we 
consider the number of points {xp,yp) S Z^\£ which are at integral distance to an embedded triangle 
£ — {{'X-a,Va], (xb,Ub )> {'X-cVcVs, we can distinguish three different embeddings of the Heronian triangle 
Ai = ( 25 , 20 , 1 5 ) . The embedding £ i = { ( , ) , ( , 25 ) , ( 1 2 , 1 6 ) } of A , y ields 1 2 points ( xp , y p ) at integral 
distance to the corners of Ai given by£i. For the embedding £2 = {(0, 0), (15, 20), (0, 20)} we obtain 1 6 
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such points, and for the embedding £3 = {(0, 0), (7, 24), (16, 12)} we obtain only 5 such points. Determin- 
ing the possible embeddings of a given Heronian triangle A = (a, b, c) is a rather easy task. W.l.o.g. we 
assume a — max{a, b, c} and xb = = ye - Since the point (xct) c) is at distance a to the point (xb ,yB), 
we have to solve the Diophantine equation 

x^ + = 

in integers. This is a well known problem. One might even store for each small number (e.g. a ^ 1 000) 
Q € N a list of the corresponding solutions (xc,yc)- Now the coordinates of the remaining point A are 
given as solutions of the equation system 

(xb -xa)^ + (Ub -Ua)^ = 



(5) 



which can be easily solved. As an algorithm for the embedding of an Heronian triangle in Z we obtain: 

Algorithm 3.4. (Embedding of an Heronian Triangle) 
input: Heronian Triangle A = (a, b, c) 
output: complete list of different embeddings of A in 
begin 

Xb = O.ijb = 

loop over the integer solutions [xcVc ] ofx^ + -y^ _ ^2 

loop over the integer solutions (xa.Ua) of equation system Q do 
output {{x A, V a), {xb,Vb], [xcVc]} 

end 

The next step is to determine the points (xp ,yp) e which are at integral distance to a given embedded 
triangle {(xA,yA), (xB,yB), {xc,Vc]}- 

Algorithm 3.5. (Enlargement of an embedded triangle) 
input: Embedded triangle £ — {(xA,yA], (x-BiUb), [xcVc]] C Z'^ 
output: complete list of points (xpiijp) G Z'^\£ which are at integral distance to £ 
begin 

loop over the integer solutions (xp,yp) of equation system (|7| do 
if (xp,'yp) ^ £ then 
output (xp,'yp) 

end 

We remark that the previous algorithms have to be implemented using an arithmetic which is able to do 
integer calculations with unlimited precision, since the occurring numbers can increase very quickly. We 
have utilized the software package CLN fS] for this purpose. 



Now we utilize the set of points given by Algorithm 3.5 to build up a graph S(£)- The vertices are given 
by the possible points (xp,yp). Two points (xp , , y p , ) and (xp^ , y P2 ) are connected by an edge if and only 
if ^(xp, — xpj )^ + (yp, — ypj)^ is a positive integer. A complete subgraph of S(£) is called a clique. A 
clique Ci is called maximal if it is not properly contained in another clique C2 of S(£)- Clearly the cliques 
of S(£) are in bijection to integral point sets T C Z^ containing £ as a subset. The same statement holds for 
maximal cliques of S{£) and maximal integral point sets T C Z^ containing £ as a subset. Thus we can use 
a clique-search package as Cliquer 121] to exhaustively generate maximal integral point sets M over Z^. 

Let us consider an example. If we apply our algorithm on the embedded triangle 

£2 ={(0,0), (15,20), (0,20)} 
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with diameter 25, we obtain a set 

{(0,28), (0,40), (0,56), (0, 132), (0,-92), (0,-16), (0, 12), (-15,20), 

(15,0), (-21,20), (105,-36), (21,20), (-48,20), (48,20), (-99,20)} 

of 16 possible points to enlarge the integral point set £2. The clique-search program Cliquer determines 
five maximal cliques which correspond to the following five maximal integral point sets: 

M, = {(0,0), (15,20), (0,20), (15,0)}, 

M2 = {(0,0), (15,20), (0,20), (0,-92), (105,-36)}, 

M3 = {(0,0), (15,20), (0,20), (0,40), (0,56), (0,-16), (-15,20), (-48,20), (48,20)}, 
M4 = {(0,0), (15,20), (0,20), (0,40), (-15,20), (-21,20), (21,20), (-48,20), (48,20), 

(-99,20), (99,20)}, and 
Ms = {(0,0), (15,20), (0,20), (0,28), (0,40), (0,56), (0, 132), (0, -92), (0,-16), (0, 12), 

(-15,20)}. 

It is interesting to have a look at the cardinalities and diameters of these maximal integral point sets. We 
have |Mi| = 4, diam(Mi ) = 25, IM2I = 5, diam(M2) = 119, IM3I = 9, diam(M3) = 96, IM4I = 11, 
diam(Mi) = 198, IM5I = 11, and diam(M5) = 224. Although we start with a point set £2 of small 
diameter, the resulting maximal integral point sets Mi may have a large diameter We are not aware of a 
formula to bound diam(M) with respect to diam(£). A second somewhat disappointing fact of our algorithm 
is, that each subset £ ' of three non-collinear points of an maximal integral point set M produces M. Thus 
our algorithm produces many identical copies of maximal integral point sets with large cardinality. We will 
deal with this fact and the isomorphism problem in the next section. 

The algorithms described in this section focus on the maximality of the integral point sets. They should 
not be used to exhaustively generate all maximal integral point sets up to a given diameter. To perform 
this task the algorithms to exhaustively generate integral point sets with or without additional properties are 
better suited, see |16, 20|, and ignore the maximality condition in the first run. All integral point sets with 
required cardinalities and small diameters can then be tested if they are maximal. 

4. Normal forms and automorphisms for integral point sets over 

In this section we aim to consider isomorphisms which preserve certain properties of maximal integral point 
sets. Since a main property of an integral point set is the set of distances between its points we only consider 
distance-preserving isomorphisms, so called isometries. In the Euclidean plane the isometries are given by 

, • -r fA fA fA ■ fA fcosd -sine\ fx 

compositions of translations Lv M + , rotations Ro : n-^ . ^ ,^ • 

' \vj \yj \vj Vsi"9 J \y 

and reflections at one of the two axes. Each isometry can be written as It,o : x i-^ t + O • x, where t e is 
a translation vector and O G M^^^ an orthogonal matrix. Next we restrict ourselves to mappings which map 
integral coordinates onto integral coordinates. Thus we have t e and O G Z^^^. Each such isometry 
It,o maps integral point sets onto integral point sets. It is easy to figure out that there are only 8 orthogonal 
matrices in Z^^^. So we define 

fr . ^2 ^ f /±i \ /±i o\ fo ±A / ±r 

Aut := o : t e Z-^, O e 



±iy ' \^ TV ' V±i y ' 

as the automorphism group of plane integral point sets. 

We call two integral point sets T and 3" isomorphic, if there exists a mapping It,o G Aut such that 
It,o (3^) = J". So our aim is to develop an algorithm which can check whether two given integral point sets 
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are isomorphic. For this purpose we want to use the technique of normal forms of discrete objects. This 
means that we have a function t which fulfills the following: If is the orbit of an integral point set CP under 
the group Aut then t( = t( J") for each J" G 0. Additionally for two integral point sets of different orbits 
the function t should have different images. Having such a function t at hand we can easily decide whether 

two integral point sets CP and CP' are isomorphic, by checking whether t(CP) — t(CP') or not. 
In order to describe such a function t we need to define a total ordering -< on Z'^: 

(1) if |a| < |c|, then we set ( ) ^ ^ 



(2) if a > 0, then we set ( ) -< 



(3) if |b| < |d|, then we set ( ) < ( 1 , and 



(4) if b > 0, then we set ( °, ) ^ 



for all Q, b, c, d G Z. We set = if and only if we have a = c and b = d. By xi ^ xz we mean 
Xi -< X2 or x^ =X2. One of the properties of this total ordering ^ is, that we have ^ x for all x e Z-^, 
so ( ^ ) ^ X is the smallest element in Z^. Using -< we can bijectively identify an integral point set CP with 



a list '^C?) of its points, which is sorted in ascending order with respect to ^. Now we extend our total 
ordering ^ onto such lists by utiUzing the lexicographic ordering. This allows us to define our normalization 
function by 

tI?) =min{£(c7(CP)) : ere Aut}. 

To obtain a finite algorithm for the determination of t(CP) we use the fact, that for every point set CP ^ the 

^0^ 



minimum list-representation £(cr(CP)) starts with y^j : 

Algorithm 4.1. (Normalization of an integral point set) 
input: integral point setV = , . . . , Prt} 
output: minimum list representation t(;P) 
begin 

champion = -C(CP) 

M,=ri ?YM,.r! °.YM,.r-' ?\m,.^-' ° 



1/ ' \0 -1/ ^ \0 1/ ' \Q -1 

for i = 1 , . . . , n do 
for j = 1 , . . . , 8 do 

tmp = £(Mj -{pi -pi,...,pn-Pi}) 
if tmp -< champion then 
champion = tmp 
return champion 

end 
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We remark that Algorithm 4.1 runs in O (n^) time. As an example we consider the two integral point 



sets from Figure [T] Their normal forms or minimum list representations are given by 



and 



-12 



-12 
-16 



-24 




-24 



respectively. 

For a given integral point set T there may exist rotation matrices M e M^^^, such that MfT) has integral 
coordinates, which are different from the eight orthogonal matrices in Z^^^. But for these matrices there is 
no guarantee for a proper extension £ D T, which is also an integral point set over Z^, such that M(£) has 
integral coordinates. Examples are given by the sets £], Ez, £3 in Section[3] This means that for a given 
maximal integral point set M over Z^ there can exist an orthogonal matrix M G M^^^, such that M(M) is 
also an integral point set over Z^, but which is not maximal. 

We may call a maximal integral point set M over Z^ strongly maximal, if such a matrix M does not exist. 
To check whether a given integral point set 7 is strongly maximal, we only have to consider all possible 
embeddings of T in Z^, which are finitely many. Another possibility is to slightly alter Algorithm 3.5 by 



which are finitely many. Another possibility is to slightly alter Algorithm ^ 
looping over the rational (instead of integral) solutions (xp,yp) of equation system ([T]i. Now the algorithm 
leads to point sets with integral distances and rational coordinates. But due to Conjecture |3.2| (which is a 
theorem for dimension m = 2), there exist embeddings with integral coordinates. 

To clear the situation with integral and rational coordinates we will have to give some facts from the 
general theory of integral point sets (without integral coordinates). So, let T be a set of points in the m- 
dimensional Euclidean space E"^ with pairwise integral distances. By 8 C T we denote an integral simplex, 
which is a set of m + 1 points, and by raboxvolmfS) we denote the m-dimensional volume spanned by 
the m + 1 points. Since the pairwise distances are integral we can write vol^(§) = q • k with q G Q and 
k being a square free integer. If voItt^(S) ^ the square free integer k is unique and we set char(S) = k, 
which we call the characteristic of §. Using this notation we can cite two results from ifTTll : 

Theorem 4.2. In an m-dimensional integral point set f all simplices § — {vq , V] , . . . , v^j with vo/m (§) 7^ 
have the same characteristic char{§) — k. 

So we can speak of the characteristic charf?) of an integral point set CP. 

Lemma 4.3. An integral m.- dimensional simplex S — {vq , , . . . , v(^} with distance matrix D — (dij ) G N 
for ^ i, j ^ m. and voZm(S) 7^ can be transformed via an isometry into the coordinates 

vo = (0,0 0), 

VI = (qi,iv/l^,0,0...,0), 

V2 = (qi.iV^. q2,2yi<2.o,...,o), 



qm.l V 1^1 > C|m,2 V 1^2, • ■ ■ , qm, 



where ki is the squarefree part of — 



voli{v'o,v\ ,...,v[)^ 



I (v^,v; ,...,v(_, ) 



qi,j G Q, and qj,j,kj ^ 0. 



MAXIMAL INTEGRAL POINT SETS OVER 1? 



9 



We remark that we always have ki = 1 . The connection between the kt and the characteristic char(T) — 
k is given by 

m 

charfT) — char(§) = k = square free part of Y\ ^i- 

i=^ 

Thus plane integral point sets T with rational coordinates are exactly those with characteristic charfT) — 



1. Due to Conjecture 3.2 plane integral point sets over correspond to plane integral point sets with 
characteristic 1. So in principle there is no need to care about the coordinates - this can still be done 
afterwards. 

There is one further transformation that maps integral point sets over onto integral point sets over 
Z^: scaling by an integral factor A. One handicap of this mapping is that the inverse mapping may lead to 
non-integral point sets. Another shortcoming is that maximal integral point sets may be mapped onto non- 
maximal integral point sets. An example is given by the maximal integral point set T = { (^o) ' (^o) ' (^4 

0\ f6\ fO\ {6 



If we scale it by a factor of 2 we obtain 2 ■ V — <^ y^j , y^j , y^j , y^j j an integral point set over 
Z^ which can be extended by the point . In contrast to this example the integral point set 3 • T — 

1 ' ' G'^) ' } might conjecture that for every maximal integral point set 

M there exists an integer A > 1 such that A • M is also maximal. 

5. Maximal integral point sets with given cardinality and minimum diameter 

From the combinatorial point of view a natural question is to ask for the minimum possible diameter 
dM(l<^, TTt) of a maximal integral point set M C Z™^ of cardinality k. If such a point set does not exist 
we set dM{k, m) — oo. Utilizing the exhaustive algorithm described in Section [5] we have obtained the 
results given in Table [T] 



k 


dM(k,2) 


corresponding point set 


4 


5 


{(0,0), (3,4), (0,4), (3,0)} 


5 


8 


{(0,0), (3,4), (0,4), (0,8), (-3,4)} 


6 


25 


{(0,0), (12, 16), (12, 9), (-12, 9), (-12, 16), (0,25)} 


7 


30 


{(0,0), (6,8), (0,8), (0, 16), (-6,8), (-15,8), (15,8)} 


8 


65 


{(0, 0), (15,36), (0, 16), (15,-20), (48, -20), (48,36), (63,0), (63, 16)} 


9 


96 


{(0, 0), (15,20), (0,20), (0,40), (0,56), (0,-16), (-15,20), (-48,20), (48,20)} 


10 


s$ 600 


{(0,0), (22, 120), (0, 120), (-27, 120), (160, 120), (182, 0), (182, 120), 

(-209, 120), (209, 120), (391 , 120)} 


11 


70 


{(0,0), (5, 12), (0, 12), (0,24), (-5, 12), (-9, 12), (9, 12), (-16, 12), (16, 12), 

(-35, 12), (35, 12)} 






{(0,0), (35, 120), (35,84), (-64,-48), (0,204), (-189,-48), (-64,252), 


12 


sc 325 


(-253,0), (-189,252), (-288,84), (-288, 120), (-253,204)} 


13 


^2046 


{(0,0), (48, 64), (0,64), (0,128), (-48, 64), (-120, 64), (120, 64), (-252, 64), 
(252, 64), (-51 0,64), (51 0,64), (-1023, 64), (1023, 64)} 



Table 1 . Minimum possible diameters of maximal plane integral point sets with given cardinality. 
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Clearly we have dM (1 , 2) — cIm (2, 2) = oo since a line I through two different points Pi and P2 with 
integral coordinates and integral distance Pi P2 contains an infinite integral point set J" = {Pi +A - (P2 — Pi ) : 
A e Z] as a subset. So the next value to determine is djvi (3, 2). Whether djvi (3,2) is finite had been an open 
question of |3 |, which was answered in |12| by determining dM(3,2) — 2066, - a diameter out of reach 
for our general exhaustive algorithm described in Section [3] But it can be easily adapted for this purpose. 
We alter Algorithm 3.1 by omitting right-angled triangles, since these obviously are not maximal. Then we 
skip Algorithm |3.4| and directly run the version of Algorithm |3.5| where we search for rational instead of 
integral solutions (xp,yp) of equation system ([T]l. If we have found the first solution (xp,i)p) for a given 
triangle A we can immediately stop our investigations on A since it cannot be a maximal integral triangle. 
Using these reductions and skipping the time consuming clique search we were able to exhaustively search 
for (strongly) maximal integral triangles over with diameter at most 15000 llT2HT8l . There are exactly 
126 such examples. Here we list the first, with respect to their diameter, ten examples, where we give the 
edge lengths and the coordinates in minimal list representation, which is unique in these cases; 



{2066,1803,505} [(0,0)'^ 

{2549,2307,1492} [(0,0)'^ 

{3796,2787,2165} [(0,0)'^ 

{4083,2425,1706} [(0,0)'^ 

{4426,2807,1745} [(0,0)'^ 

{4801,2593,2210} [(0,0)'^ 

{4920,4177,985} [(0,0)'^ 

{5044,4443,2045} [(0,0)'^ 

{5045,4803,244} [(0,0)'^ 

{5186,5163,745} [(0,0)'^ 



-336, -377)^^, (384,-2030)"^] 
-700,-2451)'^, (1100,-1008)'^] 
-387,-2760)^^, (1680,-3404)^^] 
-410,-1656)^^, (1273,2064)'^] 
-280, -2793 )'^, (376, -441 0)'^] 
-1488,-1634)'^, (1632,2015)'^] 
-473,-864)'^, (4015,1152)'^] 
-1204,-1653)'^, (2156,-4560)'^] 
-44,-240)'^, (240,4797)'^] 
-407,-624)'^, (4030,-3264)'^] 



6. Constructions for maximal integral point sets over 1? 

In this section we want to describe constructions for maximal integral point sets M of a given cardinality 
or a given shape. In most cases our constructions do not lead to integral point sets which are maximal in 
every case, but which yield candidates, which are very likely to be maximal (from an empiric point of view). 
W.l.o.g. we can assume that the origin (0, 0)^ is always contained in M. Every further point (a, b)^ meets 
+ — c^. In this case we call (a, b) a Pythagorean pair or (a, b, c) a Pythagorean triple. If additionally 
gcd(Q, b) — gcd(Q, b, c) = 1 we speak of primitive pairs or triples. Given only one Pythagorean pair (a, b) 
we can perform the following two constructions for integral point sets over 1?: 

Construction 6.1. /f (a,b) is a Pythagorean pair, then CPi (a,b) 
integral point set of cardinality 4. 

Construction 6.2. //(a, b) is a Pythagorean pair, then J'2(<i, b) :— 
is an integral point set of cardinality 5. 




We call Construction 6.1 the rectangle construction of (a,b) and Construction 6.2 the rhombus con- 
struction of (a, b). If we choose (a,b) with 2|a, 2|b then clearly (Q,b) cannot be maximal. On the 
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Other side CPi (9, 12] is a maximal integral point set although gcd(9, 12) =3. Empirically, we have ob- 
served that for primitive pairs (a,b) the point set Ti (a,b) is maximal in many, but not all cases, see 
e.g. the non maximal integral point set CPi (7,24), which can be extended to the maximal integral point 



set 



12 



For 



(3,4) both constructions J'i(a,b) and 



J'2{ti,b) yield maximal integral point sets. Empirically Construction 6.2 is a bit weaker, since it often 
happens that Ti (a,b) is maximal but J'2{ti,b) is not, as for example for (a,b) = (5, 12). For the other 
direction we have no example. We would like to mention that 72(5, 12) can be extended to the very in- 



teresting maximal integral point set M — 

>, which has an intriguing geometrical structure, see Figure 




-16^ 

. , 



'35 
. 



-35 




Figure 2. Extension of 92(5, 12) to a crab of cardinality 11. 



Definition 6.3. For positive integers a, bi , . . . , b^ we call the point set 

/ \ /±bi\ (±h^ 
^Oj\±a)'[ )-■■'[ 

a crab of order k. 



crab(a,bi , . . . ,bic) := 



We remark that the cardinality of crab(Q, bi , . . . , b^) is given by 2k + 3 and that the point set is sym- 
metric w.rt. the two coordinate axes. This point set is indeed integral if the pairs (a, bi ),..., (a, b^) are 
Pythagorean pairs. So it is very easy to construct crabs, either directly or by extending CP2 ( a, b), see Subsec- 
tion 6. 1 Empirically the extension points of J'2 ( 1, b ) very often lie on one of the two axis. An example that 
this must not be the case in general is given by the primitive pair (1480, 969), where 72(1480, 969) can be 
^0\ /1480\ /-1480\ f \ f \ /1040\ /1040\ /-1040\ /-1040 
,0j'l0 ) '[ ) ' I 969 ) ' 1-969 i ' I 462 i ' 1-462 / ' I 462 i ' I -462 



extended to 



6.1. Construction of crabs. Since many maximal integral point sets over are crabs we are interested 
in a method to construct them directly. From the general theory of integral point sets we know that integral 
point sets T over with minimum diameter consist of point sets with n — 1 coUinear points, see Figure 
[3] for 9 ^ n ^ 122 points, see IIT6I (201 . For these point sets there is an interesting connection between the 
points of the point set V and divisors of a certain number D, see | 
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The decomposition number D of an integral triangle with side lengths a, b, and c is given 

P _ (g + b + c)(a + b - c)(a - b + c) (-g + b + c) 
gcd(b2-c2 + a2,2g)2 




gs 0-2 gi go gi' CL'2 



Figure 3. Plane integral point set T with n — 1 points on a line. 
Lemma 6.5. (Decomposition lemma) 

The distances of a plane integral point set T consisting ofn points where a subset o/n— 1 points is collinear 
correspond to decompositions of the decomposition number D of the largest triangle ofT into two factors. 

Proof. We use the notation of Figure [3] and set 

1 i 

Ci = q + ^ gj for < i < s, c{ = q ' + ^ g| for ^ i < t. 
j=i j=i 

Pythagoras' Theorem yields c?_|_^ +h,^ = b?^, and c? + h.^ — b? for ^ i < s. We subtract these equations 
from each other and get 

i 

li+i + 2gi+i 2gt+i q = h^+^ - . 

j = i 

Because the g^ and the bi are positive integers we have Zg^+iq e N for ^ i < s and therefore 

2gcd{gi , g2, . . . , gs)q G N. From q + q' = go G N we conclude 2gcd(gi ,a2 gs)q' G N. With an 

analogous conclusion for the c( and g = 2gcd(gi , . . . , gg, g^ , . . . , g() we get gq e N and gq' G N. 
A last use of Pythagoras' Theorem yields for 1 ^ i ^ s and for 1 ^ j ^ t the factorization of g^K^ into a 
product of two positive integers, 

g^h^ = (gbi + gCi)(gbi - gci) = (gbj + gc-){gb- - go-). 

So we can obtain the possible values for Ct and c( by decomposing g^h.^ into two factors. 

If we are given the three side lengths g, b, and c of an integral triangle and want to determine the points 
on the side of length g so that the resulting point set is integral then we can associate b with bj, c with b(, 
and a with Y.i=^ gi + go + Y.i=^ '^i- With this we have 



12 

Definition 6.4. 

by 
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Because g can also be defined as the smallest integer with gCj G N we receive 

2a 

^ ~ gcd(b2 -c2 + a2,2a) ' 

Due to the Heron formula 1 6A^ = (a + b + c)(a + b — c)(a — b + c](— a + b + c) and the formula for 
the area of a triangle 2Aa — ah. we finally get 

2, 2 g^(a + b + c)(a + b-c)(a-b + c)(-a + b + c) 



4a2 

{a + b + c)(a + b-c)(a-b + c)(-a + b + c) 



gcd(b2-c2 + a2,2a)2 



= D. 



□ 



If we choose g = 1 and h, e N we can directly apply Lemma 6.5 to construct crabs. Let us look at 
an example. We choose g — 1 and h. = 2 • 3 • 5 = 30. The divisors of D — g^h.^ — 900 are given 
by {1,2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 25, 30, 36, 45, 50, 60, 75, 90,100, 150, 180, 225, 300, 450, 900}. If we 
have D = fl • fa, then b = ^' and c — ^' . Thus we must have f] > and fi = fi 
(mod 2) to determine the values b^ of a corresponding crab. Here we have b] = ^"^^^ — 16, ba — 
^ 40, b3 = = 72, b4 = ^^^^ = 224, and a = h = 30. This yields the integral 

point set crab(30, 1 6, 40, 72, 224) of cardinality 2-4 + 3 = 11 and diameter max |2bi^, 2a, y^b? + a^j = 

2-max{bi, a} = 448. Given the prime factorization h. = Oi=i pT^ difficult to determine the k-value 

of the resulting crab. Let us fix pi —2 and set cci = max (ai — 1,0). With this we can state 

(2ai +1).n[=2(2^i + 1) -\ 

Using h. = p'', where p is an arbitrary odd prime, we are able to produce a crab of order k for each k ^ 1 . 
Thus we have constructions for integral point sets of cardinality 2k + 3 for each k S N. To obtain small 
point sets with many points we should clearly choose integers with many divisors for h. instead. As for all 
of our constructions the maximality of the resulting integral point set is not guaranteed, but very likely. 

Construction 6.6. For a given integer h. there exists an integral point set decompose[h.) which is a crab of 
order k, where k is given by Equation 

If h, > 4 then the diameter of decompose(h,) is given by h — 1 if h. is odd and given by ^ - 2 if h IS 
even. 

Conjecture 6.7. For each integer h, the plane integral point set T = decompose[\\] is maximal if\y\ ^ 7. 

Also, the recognition of a crab is a very easy task. Given an integral point set !P over one can easily 
check whether a subset iL C T of n — 2 points is collinear by using: 



Lemma 6.8. Three points (xi ,'y i ), (x2,t)2). i^nd (x3,'y3) in are collinear if and only if we have 

0. 



xi yi 1 

X3 V3 1 



Additionally the lines through H and JXiL are perpendicular. If the point set is symmetric to these two 
lines then T is a crab. 
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Crabs are very dominating examples of maximal integral point sets over Z^. For the study of maximal 
integral point sets over 1? one might try to generalize the construction of a crab. Let us remark in this 
context, that the existence of an integral point set with coordinates 




where x, ij , z e Z is equivalent to a famous open problem, the existence of a perfect box, see Q D18]. 



So far we have only used g = 1 in Le mma 6.5 Now we want to have a look at the case g > 1 . So given 
integers g • h, and g we can apply Lemma 6.5 For two factors f i > fa with f i • fa = g^h,^ we have 



f 1 + f2 . f 1 - f2 

gbi = — - — and gct = — - — . 



The values gbi and gct are integers if and only if we have f i = fa (mod 2). Since not only the gbt's but 
also the bt's must be integers we have to require f i + fa = (mod g). Let us have an example. We choose 
gh = 672 = 2^ • 3 • 7 and g = 5. Now we look at the divisors of g^h^ = 451584 = 2^° • 3^ • 7^ and 
determine the suitable pairs (f i , fa) fulfilling 



f, -fa^g^h,^, fi>fa, f i = fa (mod 2), and fi + fa = (mod 5), 



I (784, 576), (896, 504), (11 76, 384], (1344, 336], (1536, 294), (1764, 256], (201 6, 224], 

(2304, 196), (3136, 144), (3584, 126), (4704,96), (5376, 84), (7056, 64], (8064,56), (12544,36), 

( 1 881 6, 24) , (28224, 1 6) , (32256, 14) , (75264, 6) , ( 1 1 2896, 4] } . 

The corresponding values bi are given by 

{1 36, 140, 1 56, 1 68, 1 83, 202, 224, 250, 328, 371 , 480, 546, 71 2, 81 2, 1 258, 1 884, 2824, 3227, 7527, 1 1 290} 
and the corresponding values gci are given by 

e = 1 1 04, 1 96, 396, 504, 621 , 754, 896, 1 054, 1496, 1 729, 2304, 

2646, 3496, 4004, 6254, 9396, 1 41 04, 1 61 21 , 37629, 56446| . 

Clearly the Ci cannot be integers unless h, is an integer. So let us consider the points on the left of the base 
point F of the height h. They correspond to values gci which all fulfill gct = ra (mod g), for a fixed 
ra € {1 , . . . , g — 1 }. The points on the right hand side of F correspond to the values gci fulfilling gci = — m 
(mod g). So let us choose ra = 1 . Since all elements of our candidate set 6 are congruent to ±1 modulo 5 
we obtain an integral point set of cardinality |C| + 1 = 21 : 



-9396 
5 


2304 
5 




-196 
5 




-16121 
5 




-396 
5 





-621 
5 




-896 
5 




-56446 
5 




104 
5 




504 
5 




-1496 
5 


754 
5 




-2646 
5 




1054 
5 




-3496 
5 




1729 
5 




4004 
5 




6254 
5 




14104 
5 




37629 
5 




)■} 
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After a suitable transformation and applying Algorithm 4. 1 we obtain the minimum coordinate represen- 
tation 



\ /-40\ / 64 \ /-88\ /n2\ / 144 
-168) ' I 30 j ' 1-48) ' I 66 j ' 1-84] '1-108 



-196 
147 



-920 
690 



224 \ /-288 
-168j ' 216 

1584 \ /-2176 
-1188/ ' I 1632 



320 \ / 504 \ /-560 
-240/ ' V-378y ' 420 

2660 \ /-5940 
-1995/ ' I 4455 



180 

-135 

640 
-480 ) ' 

9112 \ 

-6834 ) 



We call point sets arising from Lemma 6.5 where g > 1 and h ^ N semi-crabs, see Figure |4] for a 
drawing of our example. 



V 



Figure 4. A semi-crab of cardinality 21 and diameter 18815. 



Definition 6.9. For positive integers g, h and non-zero integers bj , . . . , we call the point set 



semi-crab (g, h, b] , . . . , bj^) := 



g 




a semi-crab of order k. 



Construction 6.10. For given positive integers g and gh, where h, ^ N, there exists an integral point set 
decompose{g, gh.) which is isomorphic to a semi-crab. 
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Conjecture 6.11. For each pair of integers gh, g the plane integral point set T — decompose{gh, g) is 
maximal if IT] ^ 7. 



Since Construction 6.6 can only produce integral point sets of odd cardinality Construction 6.10 



convenient completion. It is not difficult to implement an algorithm that exhaustively generates crabs and 



semi-crabs up to a given diameter by utilizing Lemma 6.5 



Let us have a look at the possible values for g > 1 . If we would choose g = 2 then due to 2 | g^h.^ 
all divisors of g^h.^ would be odd and we would have m = 1 (mod 2), which is not possible. Thus 2] g. 
For g = 3 the only possibility for f i + f2 = (mod 3) is f i =1 (mod 3) and f 2 = 2 (mod 3), which is 
not possible since 1-2 = 2 (mod 3) is not a square in Z3. Thus g = 5 is the first valid possibility. More 
generally we can state that if g is a prime then we have g = 1 (mod 4), since —1 has to be a square in Zg. 



6.2. Construction of integral point sets on circles. In Addition to the construction of crabs there is another 
useful construction of integral point sets of with large cardinality, see HOJ for a similar construction over 



the ring Z 



2 



Let pj = 1 (mod 4) be distinct primes over N. We consider the ring Z[i], where 



every integer pj has a unique prime factorization pj = cuj • cuj, where c denotes the complex conjugate of 
c. We may write Wj — Qj + bji, with integers aj , bj . With multiplicities vj G N we set 



and for each of the t(R) divisors of R, 



Uv] 



Ylv] 



^ Uj < Vj 



we define r[2h — Yl] 
£,s for 1 sC s s$ 2t(R) by 



=1 



Vi+Uj Vj— Ui 



ri2h-i = 1- • Ti2h for 1 ^ h. ^ t(R). With this we define vertices 



£,2 



h-k 



R ' 



1 s^hsCT(R), lce{0,l}. 



We set ris —x^ +115^ with ,i|s G Z for 1 ^ s ^ 2t(R]. We have 



2vi 2vi 



Op; 



2v, 



R^ 



This yields x^ = R — y 5, which we use to calculate 



• Its - U[ 



\^l-vl-^l+vl + '^- (2xsy s - 

\(2yl - 2yl) + i • {2x^ys - 2xtV 
4(y2_y2)2^4(^^^^_^^y^)2 



2xtyt) 

1|2 



4{vi-yi){xj-xi] 
2^(xsyt -XtDs)"^- 



■4(xs-ys -xtyt) 



Thus the distance between £,s and £,t is given by |£,s 



TlsTlt 



- £,t| = ^Ixsyt -^t^sl- Since 
■ '^sXt + VsVt + 'L{xt-ys - x^yt] 
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and 



risrit = ll<^j ll^j '^j 



n 

r 

T) -ks— kt 1 r Vi+iXi— Wi — Vi— Uj+Wi ™r.-| 

i=i 

we have that the distance between £,s and £,t is integral for every 1 ^ s,t ^ 2t(R). Additionally we can 
add the center of the circle to this point set to obtain an integral point set of cardinality 2 • t(R) + 1 having 
rational coordinates. After a suitable rotation we can achieve integral coordinates. 
So let us have an example. We choose R = 5 • 13 = 65 and successively obtain 

cui=2 + i, a)2=3 + 2i, 
rii = 65i, r|2 = 65, ri3 = — 52 + 39i, Ti4=39 + 52i, 
ri5 — —60 + 25i, ris = 25 + 60i, r|7 — —56 — 33i, rjs = —33 + 56i, 

rr r rr r 91 312. ^ 91 312. 

= -65, £,2 = 65, £.3 = ^ - ^t, £.4 = -- ^ + ^\ 
595 600. ^ 595 600. ^ 2047 3696. ^ 2047 3696. 

^'-i3""i3"'' ^'^-^ly + iy'' + 

After adding the origin (0,0)^ and applying a suitable rotation and translation we obtain the maximal integral 
point set 

"^0"\ \ /^-30\ /^-30"\ /^-63"\ /^-96\ /^-96"\ /^-126"\ /^-126^ 
,0j ' 1-32) ' I 40 / ' 1-72) ' l-16j ' I 40 / ' 1-72) ' I j ' I -32 



in minimum coordinate representation. 

Construction 6.12. For a given R which has only prime factors \) fulfilling p = 1 (mod 4) there exists an 
integral point set ci>cfe(R) consisting of 2 ■ t(R) points on a circle of radius R together with its center, where 
t(R) denotes the number of divisors ofR. 

From the above it is easy to deduce that the 2t(R) points on the circle all have pairwise even distances 
and that the diameter of this point set is given by 2R. Using this we can give another construction. 

Construction 6.13. For a given R which has only prime factors p fulfilling p = 1 (mod 4) there exists 
an integral point set ci>c/e(R) consisting of 2 ■ t(R) points on a circle of radius j, where t(R) denotes the 
number of divisors ofR. 



Conjecture 6.14. The plane integral point sets given by Construction \6. 12\ and Construction \6.13\ are max- 
imal. 

We can generalize the idea of Construction |6.13| in some way. Let t be an arbitrary integer, R be a integer 
having only prime factor fulfilling p = 1 (mod 4), and IP(R) be the integral point set given by Construction 



6.12 with radius R. By J'(R, t) we denote the point set which arises from CP(R) by scaling the point set with 
a factor \, this means dividing all distances by t. Thus J'(R, t) is a point set with pairwise rational distances 
and rational coordinates. With this we can construct a graph S containing the points of 5'(R, t) as its vertices. 
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Two vertices of S are connected by an edge, if and only if the corresponding points have an integral distance 
in TfR, t). The maximal cliques C of S correspond to integral point sets TfR, t, C). 

Construction 6.15. For a given R which has only prime factors p fulfilling p = 1 (mod 4) and a given 
integer t there exist integral point sets circZe(R, t, C) consisting of points on a circle of radius ^, where 
6 is a maximal clique of the above described graph. As an abbreviation we use circle[R,t) instead of 
c/rcZe(R, t, C). 



Conjecture 6.16. For X — & Construction 6.15 gives maximal integral point sets of cardinality t(R). 



7. Maximal integral point sets over with further conditions 



k 


dM(k,2) 


construction 


k 


dM(k,2) 


construction 


3 


= 2066 


A (2066 1803 5051 


26 


< 112895 


decompose (2 


•3-7,5) 


4 


= 5 


051(3,4) =drcle(5) 


27 


s$ 2590 


decompose (2^ 




5 


= 8 


J'2(3,4) =crab(3,4) 


28 


< 203125 


cirde (5'^ • 13) 








= crab(4,3) 


29 


1798 


decompose (2 


• 3 • 5) 


6 


= 25 


circle (5^) 


30 


< 105625 


Circle (5 • 13^ 




7 


= 30 


crab (8, 6, 15) 


31 


< 211250 


circle (54 . 132' 




8 


= 65 


circle (5-13) 


32 


< 27625 


Circle • 1 3 • 


-I -!\ 


9 


= 130 


circle (5 • 13) 


33 


< 55250 


circle (5^ • 1 3 • 


17) 


10 


625 


circle (5^) 


34 


^ 142295 


decompose (2^ 


• 3 • 7 • 1 1 , 5) 


11 


= 70 


decompose (2^ • 3) 


35 


s$ 18430 


decompose (2^ 


•3) 


12 


= 325 


circle (52 • 13) 


36 


^ 40625 


circle (55 • 13) 




13 


650 


circle (52 . 13) 


37 


^ 10366 


decompose (2^* 


•3^) 


14 


15625 


circle (5^) 


38 


< 571535 


decompose (2^^ 


•3^-7,5) 


15 


^ 8190 


decompose (2^) 


39 


< 4816895 


decompose (2^ 


•3-7,5) 


16 


^ 1105 


circle (5 • 13 • 17) 


40 


< 138125 


circle (5'* • 1 3 • 




17 


= 286 


decompose (2^ • 3) 


41 


< 73726 


decompose (iJ 




18 


sc 4225 


c'irde(52 • 13^) 


42 


< 677375 


decompose (2^ 


•32.7,5) 


19 


< 8450 


circle (52 -132) 


43 


< 4573799 


decompose {l? 


•32 - 5 •7 - 11,17) 


20 


< 8125 


circle (5^* • 13) 


44 


< 66)4998 


decompose (2^ 


•32.52.7,13) 


21 


16250 


circle (5'* • 13) 


45 


< 7001315 


decompose (2^ 


•32.72,5) 


22 


53360 


decompose (2^ • 3 • 7 • 11 , 5) 
decompose (2'* • 3) 


46 


< 64833614 


decompose {l? 


.3^* .5.7. 11,17) 


23 


< 1150 


47 


sc 7198 


decompose {1? 


.3.5) 


24 


< 5525 


circle (52 • 13-17) 


48 


< 160225 


circle (5^ •13- 


17.29) 


25 


11050 


circle (52 . 13- 17) 


49 


< 320450 


circle (5^ • 1 3 • 


17.29) 








50 


< 4064255 


decompose {iJ 


•32.7,5) 



Table 2. Best known constructions for maximal integral point sets over 7L?- in arbitrary position. 
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In Table |2] we have summarized the constructions yielding the smallest diameter of a maximal integral 
point set over Z^. Some of the values divifl^, 2) could be determined exactly by an exhaustive search, but 
for most values of k we only have upper bounds (and 301 as lower bound). In some cases, denoted by 

<, we were not able to check the maximality of the constructed point sets, since their diameter was too 
large. Looking at Table [2] we observe, that the constructions of crabs (Construction 6.6 and Construction 
6. 10 1 are very dominating. The resulting point sets contain n — 2 and n — 1 collinear points out of n 
points, respectively. So it may be interesting to study maximal integral point sets over Tl? , where no three 
points are collinear. We also say, that a point set is in semi-general position, if no three points are collinear 
By djvitl^, 2) we denote the minimum possible diameter of these point sets. We can check for this further 



condition, that no three points are collinear, by applying Lemma 6.8 Using the methods and algorithms 



described in this article, we were able to obtain some exact values and some upper bounds for djvitk, 2). 
The results are summarized in Table|3] We would like to remark that we additionally have the lower bounds 
cW(k,2) > 5525forke {n,13,14,15,17}and cWi"(k,2) > 10001 fork^ 19, k^ 20,24. 

We would like to have a closer look on the smallest known examples of maximal integral point sets in 
semi-general position consisting of an odd number of points. For cardinality 5 the two smallest point sets 
with respect to the diameter are given in minimum coordinate representation by 





-78 



-20 
21 



-20 
-99 



, (-52,-39) 



and 





-80 



-45 
28 



-45 
-108 



-96 
-40 



see Figure |5] for a drawing of the first point set. Both point sets consist of four point on a circle 6 of radii 
J Q , respectively. In each case the fifth point does not he on this circle 6, but the line through 



29. 101 and^3.53 



40 



this point and the center of C is a symmetry axis of the point set. 




Figure 5 . The smallest maximal integral point set of cardinality 5 in semi-general position. 



For cardinality 7 the two smallest examples are given by 




-285 



-180 
240 



-440\ 
-384 ) ' 



and 




-855 



-540 
720 



-1320 
1152 



-700 
240 

-2100 
720 



-880 




-2640\ 
j' 



-880 
-285 



-2640 
-855 



see Figure |6] for a graphical representation of the first example. The geometric shape of the corresponding 
two point sets is similar to the case of cardinality 5. In each case 6 points are situated on a circle C of radii 
, respectively. Again we have the symmetry axis through the seventh point and the center 



and 



3-5^-37 



5^-37 

of e. 

For cardinality 9 the two smallest examples are given by 





-504 



-64 
-252 



/612 
' l255 



612 

-759 



720 
210 



720 
-714 



/836 
' ll23 



836 > 

-627; 
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k 




construction 


k 


dM(k,2) 


construction 


3 


= 2066 


A(2066, 1803, 505) 


27 


< 305218 


circle 


52-13^-17^,8) 


4 


= 5 


?! (3,4)=circle(5) 


28 


< 203125 


circle 


5^.13] 


5 


= 120 


see Figure 5 
- — . l_J 


29 


< 9311389618298531250 


circle 


5'\S) 


6 


= 25 


circle (5^) 


30 


105625 


circle 


5^.13^) 


7 


= 925 


see Figure |6] 


31 


< 232784740457463281250 


circle 


5^°,8) 


8 


= 65 


circle(5i3) 


32 


27625 


circle 


5^1317) 


9 


= 1045 


r-| 

see Figure 7 


33 


< 412343750 


circle 


5i°-132,8) 


10 


= 625 


circle (5"') 


34 


< 152587890625 


circle 


5-) 


11 


< 2434375 


circle(5i " ,8) 


35 


< 111562500 


circle 


5'^-13*,8) 


12 


= 325 


circle(5^i3) 


36 


^ 71825 


circle 


52.13^.17) 


13 


< 60859375 


circle(5i^,8) 


37 


<363 7261 56964786376953 1250 


circle 


5^^8) 


14 


15625 


circle ( 5 '= ] 


38 


< 3814697265625 


circle 


5-) 


15 


s$ 26390 


circle(5''i3^,8) 


39 


< 10314771205 


circle 


51^-132,8) 


16 


= 1105 


circle(5i3i7) 


40 


138125 


circle 


5'*1317) 


17 


<380371 09375 


circle(5i'^,8) 


41 


<227328848 10299 1485595703 1 250 


circle 


5'»°,8) 


18 


= 4225 


circle(52i3^) 


42 


< 2640625 


circle 


5«.13^) 


19 


<950927734375 


circle(5i^,8) 


43 


<5683221202574787 1 39892578 1 250 


circle 


5^^8) 


20 


= 8125 


circle (5" 13) 


44 


< 126953125 


circle 


51°. 13) 


21 


< 659750 


circle(5'5.i3^,8) 


45 


< 7630450 


circle 


5'*-13^-172,8) 


22 


< 9765625 


circle (51°) 


46 


< 2384185791015625 


circle 


5") 


23 


<595928935571106 


circle(52^,8) 


47 


<3552013251 60924 19624328613281250 


circle 


5",8) 


24 


= 5525 


circle(5^i3i7) 


48 


< 160225 


circle 


52-1317-29) 


25 


< 4462500 


circle(5''i3'',8) 


49 


< 18854062500 


circle 


5'^-13'^,8) 


26 


< 244140625 


circle(5i^) 


50 


< 17850625 


circle 





Table 3. Best known constructions for maximal integral point sets over in semi- 



general position. 



and 





-672 



-123 
164 



-123 
-836 



-816 
340 



-816 
-1012 



-960 
280 



-960 
-952 



-1323 
-336 



see Figure 17] for a graphical representation of the first example. Here in both examples all nine points are 

and g , respectively. They both can be obtained using Construction 



52-132 



situated of circles of radii 

Now we observe that the constructions based on circles. Construction |6. 12| Construction |6T3] and Con- 
struction 6.15 are very dominating in this context. The next natural step is to also forbid four points on 
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Figure 6. The smallest maximal integral point set of cardinality 7 in semi-general position. 




Figure 7 . The smallest maximal integral point set of cardinality 9 in semi-general position. 

a circle. If no three points are on a line and no four points on a circle we speak of general position. By 
d-M (k, 2) we denote the minimum possible diameter of a maximal plane integral point set in general posi- 
tion over Z^. Without the maximality condition these point sets are also known as k2 -cluster ll22l . As we 
cannot apply our most successful constructions based on crabs and circles in this case, examples are scarce. 
For the check whether four points are situated on a circle we have a well known criterion similar to Lemma 
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Lemma 7.1. Four points (xi,ij]), ('X-iiyi), (^3,1^3), (x^,^^) in M?- are situated on a circle if and only if 



holds. 



X3 
X4 



^4 




Figure 8. The smallest maximal integral point set of cardinality 4 in general position. 



Figure 9. The smallest maximal integral point set of cardinality 5 in general position. 

In TableHwe have summarized our knowledge on dM(l^, 2). For the lower bound dM{7, 2) > 599000 
we refer to lITsl . Whether cIm (7, 2) is finite (even if we drop the maximality condition) is an open problem, 
see Q |22| . If we drop the maximality condition and the condition on the integrality of the coordinates (in 
other words characteristic one), then very recently two such examples were found, see lfT4l . The smallest 
example for k = 6 is indeed the smallest integral point set of characteristic one in general position with 
cardinality 6. We would also like to give the coordinates for the second smallest examples. For cardinality 4 
we have 



for cardinality 5 we have 




and for cardinality 6 we have 

'fO\ /-135\ f420\ / 735 \ /n55\ /1290\ 
\o) ' l^-lOOsj ' l^lOOsj ' y-392) \ 6^6) ' [^]624) 
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Figure 10. The smallest maximal integral point set of cardinality 6 in general position. 



k = ITI dM(k, 2) construction 



= 2066 
= 87 

= 165 

= 1886 

> 599000 



A{2066, 1803,505 




Figure 8 

120/ 'See Figure 



-1260 
-1023 



see Figure 



-184i 
-414 



Table 4. Best known constructions for maximal integral point sets over J? in general position. 



8. Conclusion and outlook 

We have described several constructions for integral point sets over with given cardinality that fulfill some 
further properties. Although the maximality of the resulting integral point sets cannot be guaranteed so far, 
we conjecture them to be in many cases. We have described efficient algorithms for exhaustive generation of 
maximal integral point sets over Tl? and for testing the maximality of a given integral point set. Some exact 
values of minimum diameters for given cardinalities could be obtained and several values are constructed as 
upper bounds and conjectured to be the exact values. 

It remains a task to prove the maximality of point sets resulting from some of our constructions in general. 
Clearly similar problems could be considered in higher dimensions. 
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